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ABSTRACT 

We aim to understand cloud formation in substellar objects. We combined the 
non-equilibrium, stationary cloud model of Helling, Woitke & Thi (2008; seed for- 
mation, growth, evaporation, gravitational settling, element conservation) with 
the general-purpose model atmosphere code Phoenix (radiative transfer, hy- 
drostatic equilibrium, mixing length theory, chemical equilibrium) in order to 
consistently calculate cloud formation and radiative transfer with their feedback 
on convection and gas phase depletion. We calculate the complete ID model 
atmosphere structure and the chemical details of the cloud layers. The Drift- 
Phoenix models enable the first stellar atmosphere simulation that is based on 
the actual cloud formation process. The resulting (T, p) profiles differ consider- 
ably from the previous limiting PHOENIX cases Dusty and COND. A tentative 
comparison with observations demonstrates that the determination of effective 
temperatures based on simple cloud models has to be applied with care. Based 
on our new models, we suggest a mean Teg = 1800K for the L- dwarf twin-binary 
system DENIS J0205-1159 which is up to 500K hotter than suggested in the lit- 
erature. We show transition spectra for gas-giant planets which form dust clouds 
in their atmospheres and evaluate photometric fluxes for a WASP-1 type system. 



- 2 - 



Subject headings: astrochemistry — stars: atmospheres — methods: numerical 

1. Introduction 

Today's most efficient tools to interpret the observed spectra of substellar objects, i.e., 
brown dwarfs and planets, are ID static atmosphere simulations. Comparisons with obser- 
vations, based on the solution of the radiative-transfer problem which is done in great detail 
with respect to the gas-phase opacities (Tsuji 2002, 2005; Allard et al. 2001, Ackerman & 
Marley 2001, Burrows & Sharp 1999, Barman et al. 2005), ideally yield insight into the 
atmospheric structure and chemistry providing finger prints of its evolutionary state. 

Brown dwarf and planetary atmospheres have a far more complex chemistry than stellar 
objects due to the formation of clouds, which bind chemical elements and, hence, strongly 
influences the remnant gas phase inside the atmosphere. The presence of such clouds was 
previously simplified in static model atmosphere codes. We now move a significant step 
forward by kinetically treating the chemistry of cloud formation as a phase-non-equilibrium 
process in the framework of model atmosphere simulations. Our dust model has been stud- 
ied so far for a given (T, p, t>conv) structure (T - gas temperature, p - gas pressure, fconv 
- convective velocity), and we present in this letter the ffist consistent simulation of cloud 
micro-physics and atmospheric structure which allows for the ffist time to study the feedback 
of the dust formation onto the atmospheric structure. We present our ffist consistent Drift- 
Phoenix results together with a tentative comparison to the observed spectrum of DENIS 
J0205-1159 and point out possible uncertainties in present Teg determinations (Sect. H]). 
We, however, leave aside the issue of hydrodynamical cloud formation (e.g. Showman et 
al. 2006, Knutson et al. 2007, Rauscher et al. 2007) which unavoidably has to deal with 
the turbulent closure problem (Helling et al. 2004, Helling 2007). We calculate transition 
spectra for gas-giant planets and a WASP-1 type star, and evaluate the photometric fluxes 
for the 2MASS system, for the VISIR system, and for the IRAC band systems. 

2. Method 

We combined our model describing non-equilibrium dust formation (Drift; Woitke & 
Helhng 2003, 2004; Helling & Woitke 2006; Helling, Woitke, Thi 2008) with our general- 
purpose model atmosphere code (Phoenix; Hauschildt & Baron 1999). We consider a 
stationary dust formation process, where seeds form (homo-molecular Ti02-nucleation, see 
Woitke & Helling 2004, Sect. 2.2) from a highly supersaturated gas, grow to macroscopic 
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particles of /im size, gravitationally settle into deeper layers, and eventually evaporate as the 
local temperature becomes too high for thermal stability (Woitke & Helling 2003). Coagu- 
lation, a process amongst existing grains, is omitted so far since it acts on longer time and 
size scales than the dust formation processes considered here. Typical time scales to form a 
O.l/xm Si02-grain in a gas of T = 1500K, p = 10~^ gcm~^ by continuous growth and coag- 
ulation are Tgr ^ lO'^s (Fig. 3, Woitke & Helling 2003) and Tcoag ~ 10s (Eq. 370], Rossow 
1978), respectively. The settling time scale is Xg^av ~ lO^s (Eq. 20, Woitke & Helling 2003) 
for this grain crossing a pressure scale hight {Hp = lO^cm) assuming a constant drift velocity 
t^drift ~ lOcms^^. fdrift(ct) rapidly increases for p > 10~^ gcm~^, hence Tgrav decreases inward. 

A truly static atmosphere would not contain any dust (Woitke & Helling 2004), there- 
fore we include mixing by convection and overshooting by assuming an exponential decrease 
of the mass exchange frequency in the radiative zone (Eq. 9 in (Woitke & Helling 2004) with 
P = 2.2 and r™™ = 2/ (i/pfconv)), which serves to replenish the upper atmospheres and keeps 
the cycle of dust formation running. In order to keep the computing time reasonable, we 
consider the growth/evaporation of 7 solids (MgSi03[s], Mg2Si04[s], MgO[s], Si02[s], SiO[s], 
Al203[s], Ti02[s]) made of 5 different elements for which we solve 15 possibly stiff conser- 
vation equations for 256 atmospheric layers. The whole ID atmosphere problem is solved 
iteratively in Drift-Phoenix where Phoenix provides the actual {T,p, Vconv) structure 
and Drift solves the dust moment and element conservation equations in the subsonic, 
large Knudsen number case. Drift hands back {fiV, z), Vs{z), ei(i=Mg,si,o,Ai,Ti)) with /(V, z) 
the parameterised distribution function of grain volume V at atmospheric height z, Vg the 
chemical dust composition in volume fractions for all solids s involved, and ei the remain- 
ing element abundances for all involved elements i in the gas phase (Dehn 2007). Effective 
medium and Mie theory are used to calculate the opacity of these composite and chemically 
heterogeneous dust grains in Phoenix (see Helling, Woitke & Thi (2008), Table 2). 



3. Results 

Figure [1] shows that Drift-Phoenix (T, p)-profiles for late L-type brown dwarfs ba- 
sically fall between the two limiting cases of the previous Phoenix versions DuST\EI and 
(AUard et al. 2001). Dusty-Phoenix models are hotter and Cond-Phoenix mod- 
els are cooler than our self-consistent model. We notice that, amongst the (Tpfj = 1800K, 



^AU input quantities are taken from the model presented in this paper. 
^DuSTY: Dust is considered as opacity source and element sink. 
■^COND: Dust is considered as element sink only. 
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log g = 5.0) models, the Drift-Phoenix profile has the highest temperatures inside the 
convectively active region below 1 bar. For a given local temperature, the atmospheric pres- 
sure increases inside the dust forming region at T > lOOO-ft' with increasing surface gravity 
as suggested by previous models. The result is a cloud being thermally stable at higher 
temperatures, hence sitting deeper in the atmosphere. 

Grain sizes (a), material composition Vg/Vtot, number of dust particles rid and remaining 
element abundances ei determine the gas-phase opacity and are major ingredients for the 
radiative transfer calculation. Figure [2] shows for the model with (Tgfj = 1800K, log g = 5.0) 
that the upper cloud layer is mainly made of small silicate grains with only some MgO[s] and 
Fe[s] impurities which is in accordance with our previous results (Helling et al. 2006, Helling, 
Woitke & Thi 2008). However, only a few of such silicate grains populate the cloud deck 
(Fig. [3] thick gray line). The material composition changes drastically at ~ O.lbar where 
big grains form the majority of the grain population. The dust cloud formation causes a 
local element sink (Fig. [3]). The strongest element depletion coincides with the maximum 
number of particles in pressure space. High-temperature condensable elements like Al and 
Ti are depleted in a wider atmosphere range since the associated solids are thermally stable 
over a larger temperature range than those binding Mg and Si. 



4. Application 

Brown Dwarf spectra: Drift-Phoenix was used to tentatively reproduce the observed 
spectrum of the L- dwarf twin-binary system DENIS J0205-1159 observed between 0.5/im 
and 2.6 /zm (Reid et al. 2001; Leggett et al. 2001). The determination of its mean effective 
temperatures proves difficult: Leggett et al. (2001) suggested Teg =1900K and log g = 

5.5 based on Dusty-Phoenix models and Teff = 1400 . . . 1600K from structural models. 
Golimowsky et al. (2004) suggest, depending on the system's age, an effective temperature 
between 1350K and 1700K. Vrba et al. (2004) suggest T^s =1563K with an astonishing 
precision. Brandeker et al. (2006) report similar difficulties in determining stellar parameter 
of the binary OPH 162225-240515. Based on Dusty-Phoenix and Burrows models (Burrow 
et al. 2006) they suggest T^s = (2350±150)K and Teff = (2100±100)K for their components 
A an B, respectively. 

Figure H] shows our comparison of Drift-Phoenix synthetic spectra of Teg =1800K 
with log g = 5.0 / log g = 6.0 (red solid / green dash-dot line) and Tcs =1900K with 
log g = 5.5 (blue dotted) to the Denis J0205-1159 spectrum (gray) between 0.5/im and 

2.6 fim. The synthetic spectra are normalised to the observed spectrum in the K-band at 
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A = 2.1 /iin (p). The stellar parameters have been chosen to produce the overall best fit in this 
spectral range but we cannot identify a single model that fits the whole wavelength interval. 
The J-band data are best reproduced by the (Tcff =1800K, log g = 6.0) model. The H- and 
K-band data are best reproduced by the (Teg =1800K, log g = 5.0) model. Note that the 
H band is most sensitive to chances in the surface gravity of the models. The (Teg =1900K, 
log g = 5.5) model in general produces too much fiux except for the long-wavelength edge 
of the K band. Tsuji (2005) and Gushing et al. (2007) report similar challenges for their 
models in simultaneous fitting the JHK bands. All shown models reproduce very well the 
H2O absorption edge longword 1.3 fim and 1.75 fim. 

We tentatively suggest an effective temperature of Te^ =1800K. Our tentative suggestion 
for the surface gravity is log g = 5.0 since an increase in log g broadens the H band peak 
and a decrease would narrow it, both causing even stronger discrepancies to the observed 
spectrum. It appears unlikely, according to our model, that the system's mean temperature 
should be as low as T^s ^ 1600K, since the synthetic flux would drop considerably. 
Our results furthermore suggest that L- dwarfs contain a considerable amount of dust in 
their atmospheres with a gradually changing chemical dust composition (see Fig. [21). This 
explains also why our simulations yield about the same Teg range as those in Leggett et 
al. (2001) who applied the Dusty-Phoenix limiting case, where the dust was kept in 
the atmosphere as opacity source, and consequently a high atmospheric dust content was 
assumed. The actual (T, p)-proflles, however, are very different with all its implications for 
the gas-phase and cloud formation chemistry (see Fig.[T]). Compared to Teg = 1350 — 1700K 
(Golimowsky et al. 2004), our results would then suggest an age of 2 Gyr. 
So far, our synthetic spectra are the results of a detailed micro-physical treatment of dust 
formation in phase-non-equilibrium including the formation of seed particles, their growth, 
evaporation, and the effect of gravitational settling. Grain sizes and grain size distribution, 
cloud thicknesses, dust composition, total dust content etc. are results of the simulations 
and are not prescribed as it was needed in the Dusty- and Cond-Phoenix models. Our 
simulations, however, leave space for further improvement of our synthetic spectra like for 
instance regarding the treatment of alkali line proflles in the optical (Allard et al. 2007, 
Johnas et al. 2008), the formation of liquid water droplets, or the improvement of our 
mixing modelling. 

Gas-giant transition spectra: We have calculated transition spectra for a hot and a 
medium temperature gas-giant planet with Teff = 1800K and Tcs = 1400K, respectively. 



''Gushing et al. (2007) normalise to 1.25^m in the J, 1.6/zm in the H, and 2.1/zni in the K band. 
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and log g = 3.0 in the light of a WASP-1 type host star (Tcff = 6140K, log g = 4.31, 
solar metallicity; compare Collier Cameron et al. 2007). We neglect, as a first attempt, 
the irradiation of the host star and day-night effects. For example. Barman et al. (2005) 
and Burrows et al. (2006) include these effects but do not account for the influence of 
cloud formation. Figure [5] shows D x ^P'^^^^^^/^rstar^ ^.j^g transition spectra between 0.5 /zm 
and 15 /xm corrected by the surface area ratio D — (-Rpianet/-Rstar)^ (-^planet — 8 lO'^cm — 
l.lli?j,piter, i?star = 8.63 10iOcm= 1.247?^; Mpianet = 0.5Mj„piter, " surface flux). We 
add D X Fbb(Tgg'*'^'^*)/Fbb(T|g'^'') and the photometric fluxes for the planets thermal emission 
for the VISIR and IRAC bands, and the 2MASS J H Ks bands for comparison. The transit 
contrast due to the thermal emission increases with increasing wavelength but the difference 
between the Tcs = 1800K and Tog = 1400K model become constant for A > 8fim. Our 
cloud-covered gas-giant planets produce photometric contrasts which are comparable to the 
secondary eclipse measurements by Charbonneau et al. (2005), Deming et al. (2007), Snellen 
& Covino (2007). We furthermore note a good agreement with the results on HD 209458b 
by Barman et al. (2005) despite all the uncertainties in the cloud and the day-night-effect 
modelling. 



5. Summary 

The modelling of cloud formation in substellar atmospheres has a strong impact on the 
objects temperature-pressure structure which in turn determines the cloud's chemical com- 
position, the grain size distribution function, and the cloud's location inside the atmosphere. 
These are results of our consistent simulation of substellar atmospheres and detailed non- 
equilibrium dust cloud formation. We demonstrate synthetic transition spectra for gas-giant 
planets and calculate synthetic photometric fluxes based on our consistent solution of dust 
cloud formation and radiative transfer problem. A future goal is to extend our model to 
solar-system-like planets with much cooler atmosphere. 

As we have shown for the fleld object Denis J0205-1159, stellar parameter determi- 
nations based on a comparison with synthetic spectra can vary considerable. As one con- 
sequence of this, we have set out to conduct a component-based study where e.g. cloud 
compositions, dust-to-gas-ratios, and grain sizes are compared for different cloud models 
(Helling et al. 2007, 2008). 

We thank Aleks Scholz, Andrew Collier Cameron and the anonymous referee for helpful 



''http: / /phoenix. hs.uni-haniburg.de/BrownDwarfsToPlanetsl/ 
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discussions on the paper's subject. The computer support at the School of Physics and 
Astronomy St Andrews is highly acknowledged. Most of the literature search was done with 
ADS. 
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Fig. 1. — Temperature-pressure (T, p) profiles for Tcfr = 1800K, log g = 5.0 (solar abun- 
dances) as result of different dust modelling (Dirft, Dusty, Cond) in the same model at- 
mosphere code (Phoenix). We show for comparison models with (Teg = 1800K, log g = 6.0) 
and (Teff = 1900K, log g = 5.5); see also Fig. H 




Fig. 2. — Material composition of the dust cloud in volume fractions V^/Kot {Vs - volume 
fraction of solid s, Vtot - total dust volume) and the mean particle size (a) [/xm] (thick gray) 
for the (Tcff = 1800K, log g = 5.0) model with solar element abundances. 
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Fig. 3. — Metalicity \og{e^/eu) — log(ei/eH) (e? - solar abundances of element i, H - hydrogen) 
and dust number density [cm~^] (thick gray) for the (Toff = 1800K, log g = 5.0) model. 
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Fig. 4. — Drift-Phoenix synthetic spectra in comparison to the brown dwarf twin-binary 
Denis J0205-1159 spectrum from Reid et al. (2001). The synthetic spectra are normahsed 
to the observed spectrum at A = 2.1yum, they are Gaussian convolved to 3400 samphng 
points with a spectral resolution of R=600. For (T,p) structures see Fig. [H 
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Fig. 5. — Transition spectra for two gas-giant planets seen in the light of a WASP-1 type host 
star (Teff = 6140K, logg=4.31; compare Collier Cameron et al. 2007). The synthetic spectra 
are Gauss convolved to 3000 sampling points with R= 15000. Horizontal lines indicate the 
photometric bands and the symbols give photometric fluxes for the J, H, K bands (■), the 
VISIR bands (♦) PAHl, Arlll ( 1 ), SIV ( 2 ), PAH2, SiV ( 3 ), SiC, and the IRAC bands Band 
1-4 (A). 



